Applied  Energy  108  (2013)  340-3 


Transient  modeling  and  dynamic  characteristics  of  thermoelectric  cooler 

Jing-Hui  Meng  a,  Xiao-Dong  Wang  b,c’*,  Xin-Xin  Zhang  a 


CrossMark 


a  School  of  Mechanical  Engineering,  University  of  Science  and  Technology  Beijing,  Beijing  100083,  China 

b  State  Key  Laboratory  of  Alternate  Electrical  Power  System  with  Renewable  Energy  Sources,  North  China  Electric  Power  University,  Beijing  102206,  China 
c  Beijing  Key  Laboratory  of  Multiphase  Flow  and  Heat  Transfer  for  Low  Grade  Energy,  North  China  Electric  Power  University,  Beijing  102206,  China 


H 


G  H  L  I  G  H  T  S 


.  A  complete  three-dimensional  TEC  transient  model  was  proposed. 

.  The  model  couples  heat  conduction  and  electric  conduction. 

•  Dynamic  behaviors  are  studied  at  extensive  operating  conditions. 

.  Temperature-dependent  material  properties  have  strong  effect  on  dynamic  behaviors. 
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Dynamic  characteristics  are  extremely  important  for  design  and  operation  of  thermoelectric  coolers 
(TECs).  This  paper  develops  a  three-dimensional  transient  TEC  model  based  on  the  coupling  of  heat  trans¬ 
fer  and  electric  conduction  within  semiconductors.  The  model  takes  into  account  all  thermoelectric 
effects,  including  Joule  heating,  Thomson  effect,  Peltier  effect  and  Fourier’s  heat  conduction.  For  most 
of  semiconductor  materials,  Seebeck  coefficient,  electric  conductivity  and  thermal  conductivity  are 
strongly  temperature-dependent.  Therefore,  the  present  transient  model  is  used  to  compare  dynamic 
temperature  variations  at  the  cold  and  hot  ends  with  constant  and  variable  material  properties.  Small, 
medium,  and  large  applied  currents  with  various  cooling  loads  are  adopted  as  operating  conditions. 
The  results  show  that,  at  small  currents,  constant  property  model  developed  by  this  work  can  predict 
accurately  the  dynamic  characteristics,  however,  with  the  increase  in  current,  the  temperature-depen¬ 
dence  of  properties  have  more  and  more  remarkable  effect  on  the  dynamic  temperature  variations,  espe¬ 
cially  for  high  cooling  loads.  When  the  current  is  larger  than  a  specific  value,  the  heat  transferred  from  the 
hot  end  to  the  cold  end  by  Fourier’s  heat  conduction  will  exceed  the  heat  adsorbed  at  the  cold  end  by 
Peltier  effect,  thus,  the  temperatures  at  the  cold  and  hot  ends  increase  continuously,  the  TEC  cannot  reach 
the  steady-state.  This  phenomena  is  predicted  exactly  by  the  variable  property  model,  oppositely,  the 
constant  property  model  predicts  that  the  TEC  still  supply  refrigeration. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

In  recent  years,  TECs  have  attracted  more  and  more  attention 
due  to  their  features  of  no  compressors  and  refrigerants,  large 
operating  temperature  range,  easy  to  control,  reliable  operation, 
layout  flexibility,  adaptability  and  other  characteristics  [1-3],  The 
TEC  has  been  frequently  used  for  the  cooling  of  electric  devices 
such  as  CPU,  infrared  sensor,  ice-point  reference  in  thermocouple 
thermometry,  and  refrigerators  [4],  Generally,  the  TEC  includes  a 
number  of  thermoelectric  elements,  which  are  connected  electri¬ 
cally  in  series  and  thermally  in  parallel.  The  thermoelectric  ele¬ 
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ments  are  composed  of  a  pair  of  p-  and  n-type  semiconductors. 
When  electric  current  flows  across  the  thermoelectric  elements, 
the  heat  is  transported  from  the  cold  end  to  the  hot  end  caused 
by  the  Peltier  effect. 

A  large  amount  of  works  considered  the  TEC  performance 
[5-13]  and  the  TEG  performance  [14-22]  under  steady-state 
condition,  while  few  works  [4,23-29]  investigated  its  transient  or 
dynamic  characteristics.  Actually,  transient  characteristics  are  also 
extremely  important  in  practical  TEC  operation.  First,  cooling  load 
at  the  cold  end  and/or  cooling  capacity  of  heat  sink  and  ambient 
temperature  at  the  hot  end  may  be  variable  [24],  Second,  the  start 
up  and  shut  down  characteristics  of  TECs  are  also  a  major  concern 
for  cooling  purpose  [25,26].  For  example,  for  a  fixed  cooling  load  at 
the  cold  end  and  a  fixed  cooling  capacity  of  heat  sink  at  the  hot 
end,  when  a  specific  current  is  applied  to  the  TEC,  how  much  is 
the  consuming  for  the  TEC  to  reach  the  steady-state  and  what 
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are  the  corresponding  steady-state  temperatures  of  cold  and  hot 
ends.  In  addition,  as  a  current  pulse  with  a  magnitude  of  several 
times  higher  than  the  steady-state  optimum  current  is  applied  to 
a  TEC,  the  transient  lower  temperature  than  that  reachable  at  the 
steady-state  can  be  achieved  at  the  cold  end  due  to  the  delay  of 
the  thermal  diffusion  of  the  volumetric  Joule  heat,  which  is  re¬ 
ferred  to  as  the  transient  thermoelectric  effect  [27,28], 

Up  to  now,  some  transient  models  have  been  developed  to 
understand  the  dynamic  behaviors  of  TEC.  Mostly  the  previous  tran¬ 
sient  models  were  limited  to  one-dimensional  problems  [4,24,28], 
where  the  Thomson  effect  was  ignored  and  constant  material  prop¬ 
erties  were  used  because  small  temperature  difference  between  the 
cold  and  hot  ends  was  assumed.  In  addition,  as  mentioned  by  Cheng 
et  al.  [26],  in  the  existing  models,  the  p-n  element  pair  was  simply 
treated  as  a  single  bulk  material  so  that  the  difference  in  thermal 
behavior  between  the  two  semiconductor  elements  was  not  possi¬ 
ble  to  evaluate.  Hence,  Cheng  et  al.  developed  a  three-dimensional 
model  with  p-type  and  n-type  semiconductors  as  two  separate 
parts  and  the  Thomson  effect  was  taken  into  account.  However, 
the  material  properties  were  still  assumed  to  be  constants  in  the 
Cheng’s  model. 

Recently,  we  proposed  a  general,  three-dimensional  TEC  steady- 
state  model  [29],  which  is  different  from  previous  models  [5-13]  in 
which  only  heat  conduction  equation  with  Joule  heat  and/or  Thom¬ 
son  heat  as  internal  heat  sources  is  solved.  Our  model  introduced 
the  coupling  of  heat  conduction  and  electrical  conduction,  and  con¬ 
sidered  all  the  effects  occurred  in  the  TEC.  The  model  was  used  to 
figure  out  the  performance  of  TECs  with  the  temperature-depen¬ 
dent  material  properties.  The  predictions  showed  that  the  variable 
properties  and  the  heat  losses  to  the  ambient  gas  have  significant 
effects  on  the  cooling  capacity  and  the  coefficient  of  performance 
(COP)  of  the  TEC.  Three-dimensional  temperature  distributions 
within  semiconductors  was  observed  and  it  became  more  remark¬ 
able  at  large  temperature  differences  and  high  currents. 

The  purpose  of  this  work  reaches  the  following  targets:  (1) 
extending  our  previous  steady-state  model  [29]  to  a  transient 
one;  (2)  employing  the  developed  transient  model  to  investigate 
the  dynamic  characteristics  of  TECs  under  various  operating  condi¬ 
tions.  The  dynamic  behaviors  for  TECs  with  temperature-dependent 
properties  are  analyzed  and  compared  with  those  constant 
properties. 


2.  Model  development 

The  schematic  of  a  TEC  device  containing  one  semiconductor 
element  pair  (referred  to  as  TEC  element  hereafter)  is  shown  in 
Fig.  1.  A  TEC  device  consists  of  a  number  of  TEC  elements.  In  con¬ 
sideration  that  the  thermal  characteristics  among  the  TEC  ele¬ 
ments  are  periodic,  hence,  only  one  element  is  considered  here. 
The  TEC  element  includes  a  pair  of  p-  and  n-type  semiconductor 
columns,  three  metallic  connectors,  and  two  electrically  insulating 
ceramic  plates.  The  ceramic  plate,  connector,  and  semiconductor 
columns  have  thicknesses  of  H0,  H,,  and  H2,  respectively.  The  semi¬ 
conductor  columns  have  square  cross-section  with  the  side  length 
of  L2.  The  distance  between  p-  and  n-type  semiconductors  is  Lv  A 
heat  sink  with  the  thickness  of  Hs  is  attached  on  the  hot  end  of  the 
TEC  element  for  heat  dissipation.  The  geometric  parameters  of  the 
TEC  element  simulated  in  the  present  work  are  as  follows, 
Li  =  0.2  mm,  I2  =  0.5  mm,  H0  =  0.2  mm,  Hi  =  0.1  mm,  H2  =  1.0  mm, 
and  Hs  =  0.2  mm. 

When  applied  current  flows  from  the  p-type  semiconductor  to 
n-type  semiconductor,  the  holes  in  the  p-type  semiconductor  and 
the  electrons  in  the  n-type  semiconductor  migrate  from  the  cold 
end  to  the  hot  end,  the  corresponding  Peltier  heats  will  be  gener¬ 
ated  at  the  interface  between  connectors  and  semiconductors. 


The  heat  is  adsorbed  at  the  cold  end  and  liberated  at  the  hot  end 
caused  by  the  Peltier  effect,  forming  a  temperature  difference 
AT  =  TH  -  TL.  With  the  temperature  difference,  the  heat  will  be 
transferred  from  the  hot  end  to  the  cold  end  due  to  Fourier’s  heat 
conduction.  In  the  role  of  the  current  and  temperature  gradient, 
the  Joule  heat  and  the  Thomson  heat  will  be  generated  within 
semiconductors.  Thus  the  final  cooling  capacity  of  the  TEC,  Ql,  is 
determined  by  cooperative  effect  of  Peltier  heat,  Fourier’s  heat 
conduction,  Joule  heat,  and  Thomson  heat. 

2.1.  Governing  equations 


The  transient  TEC  model  is  developed  from  our  previous  steady- 
state  model  [29]  by  adding  a  transient  term.  The  basic  governing 
equations  include  the  energy  equations  and  the  electric  potential 
equations.  The  energy  equations  of  connectors,  p-  and  n-type  semi¬ 
conductors,  ceramic  plates,  and  heat  sink  are  as  follows: 

(PCp)conn  ^  =  V  ■  (AconnVT)  +  -  (ic0J  •  VT 

0) 

(PCp)p§=V.(2pVT)+£-/y.VT 

(2) 

C PCp)n §  =  V  •  (X„VT)  +L-PJ.VT 

(3) 

(PCP)cer^=V-(2CerVT) 

(4) 

(PCP)stakg  =  V-(2sinkVT) 

(5) 

where  p  is  the  density,  cp  is  the  specific  heat,  X  is  the  thermal 
conductivity,  er  is  the  electric  conductivity,  and  [i  is  the  Thomson 
coefficient.  The  subscripts  conn,  p,  n,  cer,  and  sink  denote  connector, 
p-  and  n-type  semiconductors,  ceramic  plate,  and  heat  sink,  respec¬ 
tively.  The  first  terms  on  the  right  side  in  Eqs.  (1)— (5)  denote  Fou¬ 
rier’s  heat  conduction,  the  second  and  third  terms  denote  internal 
heat  sources  due  to  Joule  heating  and  Thomson  effect,  respectively. 
J(x,  y,  z )  is  the  local  current  density,  which  is  assumed  to  be  constant 
and  equal  to  //A  in  the  previous  models  [4-13,23-28],  where  I  is  the 
total  applied  current,  and  A  is  the  cross-sectional  area  of  semicon¬ 
ductors.  The  Thomson  coefficient  can  be  derived  from  Seebeck  coef¬ 
ficient,  or: 


where  a  is  the  Seebeck  coefficient.  The  electric  potential  is  the  driv¬ 
ing  force  of  the  electrons  and  holes  in  the  semiconductor,  which  can 
be  obtained  by  solving  the  following  equation: 

V  •  (cr(V(/>  -  <xVT))  =  0  (7) 

where  <j>  is  the  electric  potential,  ocVT  is  Seebeck  electromotive  force 
coming  from  Seebeck  effect.  Once  #  is  determined,  the  electric  field 
can  be  calculated  by  the  following  equation: 

E  =  -V</>  +  aVT  (8) 

Finally,  the  current  density  vector  in  Eqs.  (1 )— (3)  can  be  calcu¬ 
lated  as  follows: 
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Fig.  1.  The  schematic  of  TEC  element. 


(9) 


2.2.  Initial  conditions  and  boundary  conditions 

The  initial  conditions  are  as  follows.  The  temperature  of  the  TEC 
element  is  uniform  and  equals  to  the  ambient  temperature  of 
T  =  300  K,  and  the  initial  electric  potential  of  the  TEC  element  is  0. 

The  adopted  boundary  conditions  are:  on  the  top  of  the  heat 
sink,  heat  is  transferred  to  the  ambient  by  convective  heat  transfer, 
or 


-hsnQ  =M^2(Ho+Hl)+«2+Hs-rj  (10) 

lz=2(H0+H,)+H2+Hs 

where  heff  is  the  equivalent  convective  heat  transfer  coefficient, 
which  will  be  determined  in  the  next  section. 

On  the  bottom  of  the  TEC  element,  a  constant  heat  flux  is  as¬ 
sumed  as  follows: 


On  the  side  surfaces  of  the  TEC  element,  the  adiabatic  boundary 
conditions  are  specified,  or: 


On  the  internal  interfaces  between  adjacent  materials,  the  tem¬ 
perature  and  heat  flux  are  assumed  to  be  continuous.  The  electric 
boundary  conditions  for  the  TEC  element  are  as  follows:  a  fixed 
current  is  applied  to  the  TEC  element  and  a  zero  electric  potential 
is  specified  on  the  side  surface  of  connector  with  x  =  2L^+  2 L2 
(Fig.  1)  or: 

f|x=0  =  constant  0  <y  <  L2,H0  <  z  <  H0  +  Hi  (13) 

4>\x-2l,+2l2  =  0  0  <  y  <  L2, Ho  <  z  <  H0  +  H,  (14) 

On  the  other  surfaces,  the  current  cannot  flow  out  of  the  TEC 
element,  thus,  we  have: 

Jn  =  0  (15) 


2.3.  Material  properties 

The  n-type  and  p-type  semiconductors  are  assumed  to  be  Bi2(- 
Teo.94Seo.o6)3  and  (Bio.25Sb0.75)Te3.  Minnich  et  al.  [3]  proposed  that 
the  properties  of  some  thermoelectric  materials  are  strongly  tem¬ 
perature-dependent,  hence  for  the  purpose  of  constructing  a  gen¬ 
eral  TEC  model,  the  variable  properties  should  be  considered.  In 
the  present  work,  the  thermal  conductivity,  electric  conductivity 
and  Seebeck  coefficient  are  all  assumed  to  be  the  functions  of  tem¬ 
perature  as  follows: 


Table  1 

Parameters  in  Eqs.  (16)-(18)  for  thermoelectric  materials  (T0  =  300  K). 


Material  Thermal  conductivity  Electric  resistivity  Seebeck  coefficient 

_ MWm-MC-1)  A,  (!<-') _ AAV.-2)  ft(Qm)  B,  (K"1)  B2  (K~2)  ^[Vr1)  C,  (K1)  C2  (K~2) 

p-Type  1.472  -1.29  x  10  3  1.35  x  10  5  8.826  x  10  6  5.88  x  10  3  8.93  x  10  6  2.207  x  10  4  1.55  x  10  3  -3.15  x  10  6 

n-Type  1.643  -9.80  x  10  4  1.56  x  10  5  8.239  x  10  6  4.70  x  10  3  2.67  x  10  6  -2.230  x  10  4  5.62  x  10  4  -4.65  x  10  15 

Copper  400  0  0  1.700  x  IQ-9  0  0  6.500  x  IQ-6  0  0 
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Table  2 

Properties  for  AIN  ceramic  plate  (T=  300  K)  [22], 


Table  3 

Properties 


of  heat  sink  system. 


Parameter 


Parameter 


Unit  Source 


p«r  3250  kgm~3 

260  wr'm4 

Cp.cer  7.6  X  105  J  kg-1  K  ' 


Kff 


Fig.  2.  The  schematic  of  heat  sink. 


Fin 

Height  Hf 

Thickness  t/ 

Length  Lsink 

Fin  to  fin  spacing  Dg 

Base  width  Wsin 

Base  height  bf 

Coolant  (air) 

Thermal  conductivity  Xa 
Density  p„ 

Viscosity  p„ 


14.0  mm  [4] 

0.2  mm  [4] 

0.5  mm 

0.2  mm  [4] 

0.2  mm  [4] 

26.3  xlO3  Wm-'K-1  [23] 

1.1614  kgm-3  [23] 

184.6  xlO-7  Nsm2  [23] 

1007  J  kg1  K-1  [23] 


Qh  =  fWKCk  -  Too)  (19) 

where  Qh  is  the  total  heat  dissipation  of  real  heat  sink,  T*°pk  is  the 
top  temperature  of  equivalent  plate  heat  sink,  and  A  is  the  base  area 
of  heat  sink,  which  equals  the  area  of  ceramic  plate,  or 
A  =  Lsink  x  Wsink  =  L2  x  (2Lt  x  2L2). 

For  a  finned  structure  heat  sink,  Mereu  et  al.  [31]  measured 
experimentally  the  relationship  between  the  heat  dissipation  with 
heat  sink  structure  and  operating  conditions,  expressed  as: 


A(T)  =  A(T0)  [1  +  Ai  (T  -  T0)  +  A2  (T  -  T0)2]  (1 6) 

^  =  P(T)  =  P(T0)[1  +Bi(T-  To)  +  B2(T  -  T0)2}  (17) 

oc(T)  =  a(T0)[l  +  C,(T  -  T0)  +C2(T-  T0f$  (18) 

where  T0  =  300  K  is  the  reference  temperature,  /.(T0),  p(T0),  and  a(T0) 

are  the  thermal  conductivity,  electric  resistivity,  and  Seebeck  coef¬ 
ficient  at  reference  temperature,  respectively:  Ai,  A2,  B i,  B2,  Ci, 
and  C2  are  the  fitted  parameters  determined  by  experimental  data 
as  listed  in  Table  1. 

The  connector  and  heat  sink  are  made  of  copper  and  their  prop¬ 
erties  are  assumed  to  be  constant  due  to  weak  thermoelectric  ef¬ 
fect.  AIN  is  selected  as  material  of  ceramic  plates  and  its 
properties  are  listed  in  Table  2  [30], 

3.  Simplification  of  heat  sink 

For  a  TEC,  it  is  desirable  to  achieve  a  high  cooling  capacity  and  a 
large  coefficient  of  performance  (COP)  through  selecting  a  proper 
applied  current  under  the  same  heat  sink  cooling  condition.  Thus, 
for  a  fair  comparison,  the  heat  sink  is  maintained  at  the  same  oper¬ 
ation  conditions,  which  means  that  the  equivalent  convective  heat 
transfer  coefficients  heff  in  Eq.  (10)  are  the  same  for  all  transient 
simulations. 

The  most  common  configuration  of  heat  sink  used  in  current 
applications  is  of  the  fin  structure  with  the  geometry  shown  in 
Fig.  2.  The  heat  sink  comprises  a  series  of  parallel  fins  with  height 
Hf,  thickness  t/,  and  length  Lsink  =  L2.  Each  fin  is  spaced  by  a  gap  Dg 
and  mounted  on  the  heat  sink  base  with  an  area  of  Lsink  x  Wsink.  Air 
is  employed  as  the  coolant.  The  geometric  parameters  of  the  heat 
sink  and  the  air  properties  are  shown  in  Table  3  [31,32],  In  order 
to  simplify  the  calculation,  the  heat  sink  with  real  structure  is  sim¬ 
plified  to  a  flat  plate  with  the  same  area,  which  was  also  adopted 
by  Chen  et  al.  [33],  In  order  to  ensure  the  same  heat  dissipation 
for  the  equivalent  and  real  heat  sinks,  the  equivalent  convective 
heat  transfer  coefficient  he ff  between  the  equivalent  heat  sink 
and  the  ambient  must  meet: 


PaWsmkHf  Dg  A p 

i  +4  «a5T 


Too) 


(20) 


where  pc,  pc,  cpc  are  the  density,  viscosity,  and  specific  heat  of  air, 
Ap  is  the  pressure  drop  across  the  heat  sink,  Tmax  is  the  maximum 
temperature  observed  in  the  heat  sink.  Copper  is  adopted  as  the 
material  of  heat  sink,  and  its  thermal  conductivity  is  400  W  m-1  K  1 
which  is  much  higher  than  that  of  ceramic  plate,  hence  the  internal 
temperature  difference  of  heat  sink  is  expected  to  be  very  small.  The 
simulations  verify  the  guess  by  the  fact  that  the  internal  tempera¬ 
ture  difference  is  less  than  0.2  K,  therefore,  we  have  T“npk  =  Tmax. 
The  equivalent  heat  transfer  coefficient  can  be  calculated  from 
Eqs.  (19)  and  (20),  or: 


Qh  Pafy  Dl  Ap 

1+5-  i2^ia  Lginka 


(21) 


In  the  present  work,  with  Ap  =  1.3  Pa,  heff  =  7712  Wm  2  K  1. 
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(b)  1  =  0.5  ;  (c) 7=  1.0 A. 


load  test  is  to  avoid  the  errors  in  measuring  the  magnitude  of  the 
cooling  load  and  the  heat  losses  with  it  which  could  lead  to  uncer¬ 
tainties  in  the  testing  conditions  [26],  Fig.  4  compares  the  measure¬ 
ment  data  regarding  the  dynamic  temperature  variations  at  the  cold 
and  hot  ends  with  numerical  predictions  by  the  present  model  at 
/  =  0.3, 0.5  and  1  A,  respectively.  It  is  found  that  the  predictions  agree 
well  with  the  experiments,  the  maximum  temperature  difference 
between  numerical  predictions  and  the  measurement  data  is  less 
than  1 .5  K  at  /  =  0.3  A,  2.0  K  at  /  =  0.5  A,  and  4.5  K  at  /  =  1  A,  respec¬ 
tively,  hence,  the  present  transient  model  can  accurately  analyze 
the  dynamic  characteristics  of  TECs. 


5.  Simulation  cases 

Thermal  conductivity,  electrical  conductivity,  and  Seebeck  coef¬ 
ficient  of  semiconductor  materials  are  all  temperature-dependent. 
Our  previous  study  [29]  has  shown  that  these  properties  signifi¬ 
cantly  affect  the  cooling  capacity  and  the  COP  of  TECs  when  a  large 
temperature  gradient  occurs  within  semiconductor  materials.  For 
the  n-type  and  p-type  semiconductors  of  Bi2(Teo.94Seo.o6)3  and 
(Bio.25Sb0.75)Te3,  with  a  small  applied  current,  there  is  no  marked 
difference  in  the  cooling  capacity  and  COP  predicted  by  constant 
property  model  and  variable  property  model,  however,  as  the  cur¬ 
rent  is  larger  than  1.2  A,  ignoring  the  property  change  overesti¬ 
mates  the  cooling  capacity  and  COP.  Especially,  note  that  the 
constant  property  model  significantly  overestimates  the  maximum 
working  current  of  TEC,  which  is  defined  as  the  corresponding 
maximum  applied  current  when  the  cooling  capacity  is  greater 
than  zero. 

In  the  present  work,  the  TEC  geometry  and  the  convective  heat 
transfer  coefficient  of  heat  sink  remain  unchanged,  with  the  same 
cooling  load  at  the  cold  end,  the  response  time  of  temperatures 
(defined  as  the  consuming  time  for  the  temperature  of  the  cold 
and  hot  ends  reaches  the  steady-state),  the  dynamic  temperature 
variations,  as  well  as  the  steady-state  temperatures  and  tempera¬ 
ture  difference  at  the  cold  and  hot  ends  are  closely  dependent  on 
the  applied  current.  At  different  applied  currents,  the  Joule  heat 
and  Thomson  heat  generated  within  semiconductors  have  large 
difference,  which  leads  to  a  remarkable  difference  in  temperature 
distributions  of  semiconductors.  Hence,  it  can  be  expected  that 
the  variable  property  model  could  predict  the  TEC  dynamic  charac¬ 
teristics  more  accurately  when  a  large  temperature  gradient  occurs 
within  semiconductors.  Therefore,  this  work  focuses  on  the  com¬ 
parisons  of  the  dynamic  characteristics  of  the  TEC  predicted  by 
constant  and  variable  property  models  under  different  applied  cur¬ 
rents  and  different  cooling  loads  at  the  cold  end.  The  parameters 
used  in  the  simulations  are  listed  in  Table  4,  where  f  =  0.6  A, 
1.5  A,  2.0  A  and  2.8  A  denote  typical  small,  medium  and  large  ap¬ 
plied  currents,  respectively,  and  qL  is  varied  from  14,285.7  W  m  2 
to  1 71,428.6  Wm  2  denoting  typical  low,  medium  and  high 
cooling  loads  at  the  cold  end. 


5.7.  At  small  applied  current 


4.  Model  validation 

The  model  used  uniformly  distributed  elements  in  the  x,  y  and  z 
directions.  The  grid  independence  was  examined  in  preliminary 
test  runs.  The  independence  of  the  time  step  size  was  also  exam¬ 
ined  with  At  =0.1,  0.5,  1,  and  2  s,  as  shown  in  Fig.  3,  and  then 
At  =  1  s  was  selected. 

Cheng  et  al.’s  experimental  data  [34]  are  used  to  validate  the 
present  transient  model.  A  zero-cooling-load  test  is  adopted  in  the 
present  comparison.  The  reason  for  performing  the  zero-cooling- 


Fig.  5  shows  the  dynamic  characteristics  predicted  by  the  con¬ 
stant  and  variable  property  models  under  different  cooling  loads 
at  small  applied  current  of  7  =  0.6  A.  It  can  be  observed  that  the  re¬ 
sponse  time,  the  temperature  difference  between  the  cold  and  hot 
ends,  and  the  temperatures  at  the  cold  and  hot  ends  predicted  by 
both  models  are  almost  identical.  These  can  be  explained  by  the 
fact  that  at  small  current  the  Joule  heat  and  the  Thomson  heat  gen¬ 
erated  inside  the  semiconductors  are  very  low,  leading  to  a  small 
temperature  gradient  within  the  semiconductors,  and  hence  a 
invisible  property  change. 
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Table  4 

Cases  simulated  in  the  present  work. 

Cases  I  (A)  qL,  (W  nr2)  qu  (W  nr2)  qu  (W  nr2)  qM  (W  m  2)  q ^  (W  m^2) 


Low  current  0.6  0  14285.7  28571.4  42857.1  71428.6 

Medium  current  1.5  0  57142.9  85714.3  114285.7  142857.1 

2.0  0  28571.4  71428.6  142857.1  171428.6 

Urge  current  2.8  0  28571.4  71428.6  142857.1  171428.6 


Fig.  5  shows  that  the  cold  end  temperature  decreases  but  the 
hot  end  temperature  increases  with  the  time  increased.  However, 
both  gradually  approach  to  the  steady-state.  With  the  increase  in 
the  cooling  load,  the  temperatures  at  the  cold  and  hot  ends  are 
both  increased,  but  the  temperature  rising  rate  of  the  cold  end  is 
much  higher  than  that  of  the  hot  end,  resulting  in  a  significant 
reduction  in  the  temperature  difference  between  the  cold  and 
hot  ends.  This  is  responsible  for  that  with  the  same  heg  and  I,  only 
a  low  temperature  difference  can  reduce  Fourier’s  thermal  conduc¬ 
tion,  hence  the  TEC  maintains  a  high  cooling  capacity.  It  is  noted 
that,  since  the  temperature  difference  decreases  and  the 
steady-state  temperature  at  the  cold  end  is  closer  to  the  initial 
temperature  of  300  K  under  high  cooling  load,  the  TEC  response 
time  is  significantly  reduced  as  the  cooling  load  increases. 


5.2.  At  medium  applied  current 

Fig.  6  shows  the  dynamic  temperature  variations  at  the  cold  and 
hot  ends  predicted  by  the  constant  and  variable  property  models  at 
/  =  1.5  A  with  various  cooling  loads.  With  zero  cooling  load,  the  dy¬ 
namic  temperature  variations  at  the  cold  end  are  almost  the  same 
with  that  at  the  hot  end,  only  a  slightly  difference  occurs  when 
t  >  200  s.  The  two  models  both  predict  the  temperature  overshoot 
phenomenon  at  the  hot  end,  for  example,  the  variable  property 
model  indicates  that  the  hot  end  temperature  increases  rapidly 
from  the  initial  stage  of  300  K,  then  decreases  slowly  after  reaching 
the  maximum  value  of  348.1  K,  finally  approaches  to  the  steady  va¬ 
lue  of  342  K  at  the  time  of  1200  s. 

The  constant  and  variable  property  models  predict  that,  at 
/=  1.5  A  the  temperatures  at  the  cold  and  hot  ends  both  increase 
with  the  increase  in  the  cooling  load,  which  has  been  observed  at 
small  current  of  0.6  A  (Fig.  5).  However,  at  medium  current  of 

I . 5  A,  the  difference  in  dynamic  temperature  variation  predicted 
by  both  models  occurs  and  the  difference  become  more  remark¬ 
able  at  high  cooling  load.  The  constant  property  model  underesti¬ 
mates  the  temperatures  at  the  cold  and  hot  ends.  In  details,  the 
steady-state  temperatures  at  the  hot  end  and  cold  end  predicted 
by  the  constant  property  model  under  qL  =  85,714.3  Wm~2  is 
347.8  K,  and  281.6  K,  which  underestimates  by  6.8  K  and  6.7  K 
compared  to  the  variable  property  model,  respectively.  Under 
qL  =  142,857.1  W  m  2,  however,  the  underestimations  increase  to 

II. 0  K  and  14.7  K  at  the  hot  end  and  cold  end,  respectively.  Fig.  7 
shows  the  temperature  distributions  within  p-  and  n-type  semi¬ 
conductors  predicted  by  both  models  under  qL  =  142,857.1  W  m  2, 
and  t  =  1000  s.  It  can  be  seen  that  the  temperature  predicted  by  the 
constant  property  model  is  lower  than  that  by  the  variable  prop¬ 
erty  model.  According  to  Eqs.  (2)  and  (3),  the  reason  can  be  ex¬ 
plained.  Because  the  electrical  conductivities  of  Bi2(Te0.94Se0.06)3 
and  (Bi0.25Sbo.75)Te3  have  negative  temperature  effect,  the  constant 
property  model  underestimates  the  Joule  heat  generated  within 
the  semiconductors,  leading  to  a  lower  temperature  distribution 
(Fig.  7). 

Note  that  with  high  cooling  loads  of  qL  =  1 14,285.7  W  m  2  and 
142,857.1  W  m  2,  a  undershoot  phenomenon  of  cold  end  temper¬ 
ature  is  observed  (Fig.  7c  and  d).  The  cold  end  temperature  de¬ 
creases  from  300  K  to  the  minimum  value  at  the  initial  stage, 


then  increases  until  reaching  the  steady-state.  This  phenomenon 
can  be  explained  as  follows:  the  Peltier  cooling  occurs  at  the  cold 
end  of  the  semiconductors,  while  the  Joule  heating  and  Thomson 
effect  occurs  uniformly  across  the  semiconductors  columns.  Dur¬ 
ing  the  dynamic  process,  the  cooling  at  the  cold  end  occurs  before 
the  Joule  heat  and  Thomson  heat  reach  this  end  and  a  temporary 
temperature  overshoot  at  the  cold  end  can  be  observed.  However, 
for  small  qL,  there  is  a  large  temperature  differences  between  the 
cold  and  hot  ends,  leading  to  a  enhanced  Fourier’s  back  heat  con¬ 
duction  from  the  hot  end  to  the  cold  end,  hence  the  undershoot 
does  not  occur. 

Fig.  8  shows  the  electric  potential  difference  and  COP  (=  QJ 
P  =  QlW  where  Ql  =  qiA  is  the  cooling  capacity  of  the  TEC,  P  is 
the  electric  power,  V  is  the  electric  potential  difference)  predicted 
by  the  constant  and  variable  property  models  at  /=  1.5  A.  It  can  be 
seen  that  due  to  the  negative  temperature  effect  of  electric  conduc¬ 
tivity  for  Bi2(Teo.94Seo.o6)3  and  (Bi0.25Sbo.75)Te3  materials,  the  vari¬ 
able  property  model  predicts  higher  electric  potential  difference 
through  the  TEC.  Therefore,  more  electric  power  is  needed  for  the 
variable  property  model,  leading  to  a  reduced  COP. 

Fig.  9  shows  the  dynamic  temperature  variations  at  the  cold  and 
hot  ends  at  /  =  2.0  A  with  qL  =  0  W  m  2  and  142,857.1  W  m  2.  It  is 
noted  that  at  the  larger  applied  current  of  1  =  2.0  A,  there  is  a 
remarkable  difference  between  the  constant  and  variable  models 
even  at  zero  cooling  load  due  to  significantly  enhanced  Joule  heat. 
With  high  cooling  load  qL  =  142857.1  W  m  2,  the  temperature  at 
the  cold  end  predicted  by  the  variable  property  model  does  not  de¬ 
creases  but  increases  continuously  to  steady-state  value  of  392.4  K, 
which  is  even  higher  than  the  temperature  at  the  hot  end  predicted 
by  constant  property  model.  The  results  again  confirm  that  tem¬ 
perature-dependent  properties  have  significant  effect  on  dynamic 
characteristics  of  the  TEC. 


constant  and  variable  property  models  at  /  =  0.6  A. 
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variable  property  models  at  I  =  1.5  A. 


5.3.  At  large  applied  current 

Fig.  10  shows  the  dynamic  temperature  variations  at  the  cold 
and  hot  ends  at  /  =  2.8  A  with  qL  =  0  W  nr2  and  142,857.1  W  nr2, 
respectively.  At  large  applied  current,  the  temperature  gradient 
within  the  semiconductors  is  significantly  elevated  due  to  the 
strong  Joule  heating,  resulting  in  strong  effect  of  material  proper¬ 
ties.  Consequently,  dynamic  temperature  variations  at  the  cold 
and  hot  ends  predicted  by  constant  property  model  remarkably 
differ  from  those  predicted  by  variable  property  model.  The  con¬ 
stant  property  model  predicts  TL  =  295.0  K  and  TH  =  430.2  K  at 


qL  =  0  W  nr2,  while  TL  =  348.3  K  and  TH  =  438.8  K  at 
qL  =  142,857.1  W  m  2  when  the  TEC  element  approaches  the  stea¬ 
dy-state.  However,  the  variable  property  model  indicates  that  the 
temperatures  at  the  cold  and  hot  ends  increase  rapidly,  and  it 
can  never  reach  the  steady-state 

Using  our  previous  steady-state  model  [29],  the  I-Qy  perfor¬ 
mance  curves  with  fixed  temperatures  at  the  cold  and  hot  ends, 
TL  =  295.0  K  and  TH  =  430.2  K,  as  well  as  TL  =  348.3  K  and 
Th  =  438.8  K,  can  be  predicted  as  shown  in  Fig.  11.  The  black  lines 
are  predicted  by  constant  property  model,  and  red  dash  dot  lines 
are  predicted  by  variable  property  model.  Note  that  the  hollow  cir¬ 
cle  in  /-Ql  curve  (constant  property  model,  TL  =  295.0  K  and 


property  models  at  7  =  1.5  A 
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Th  =  430.2  K)  denotes  the  steady-state  point  with  1  =  2.8  A  and 
QL  =  0.1  W  (gL  =  QiJA  =  142,857.1  Wm  2),  the  solid  circle  in  f-d. 
curve  (constant  property  model,  TL  =  348.3  K  and  TH  =  438.8  K) 
denotes  the  steady-state  point  with  /  =  2.8  A,  Q.  =  0  W.  The  result 
indicates  that  the  transient  model  developed  in  this  work  can 
predict  well  the  dynamic  response  characteristics  and  the  final  stea¬ 
dy-state.  The  variable  property  model  predicts  that  the  effective 


working  current  ranges  from  0.49  A  to  2.40  for  fixed  TL  =  295.0  K 
and  TH  =  430.2  K,  and  from  1.09  A  to  1.94  A  for  fixed  TL  =  295K 
and  Th  =  430.2  K.  Apparently,  for  the  transient  cases,  applied  current 
1=2.8  A  is  higher  than  the  maximum  working  current.  Under  this 
circumstance,  the  heat  transferred  from  the  hot  end  to  the  cold 
end  by  Fourier’s  heat  conduction  is  still  larger  than  the  heat  ad¬ 
sorbed  at  the  cold  end  by  Peltier  effect  even  at  zero  cooling  load, 
thus,  the  temperatures  at  the  cold  and  hot  ends  both  increase  result¬ 
ing  the  failure  of  TEC  refrigeration. 


6.  Conclusions 

This  work  develops  a  complete  three-dimensional  transient  TEC 
model,  which  introduces  the  coupling  of  the  heat  transfer  and  elec¬ 
tric  conduction  within  semiconductors,  and  the  transient  TEC  mod¬ 
el  developed  in  this  paper  has  not  been  found  in  the  existing 
literatures.  The  developed  model  is  used  to  investigate  the  dy¬ 
namic  characteristics  of  TECs.  A  miniature  TEC  element  with  a 
finned  structure  heat  sink  is  investigated  and  Bi2(Te0.94Se0.06)3 
and  (Bi0.25Sb0.75)Te3  are  selected  as  the  n-type  and  p-type  thermo¬ 
electric  materials,  respectively.  The  present  model  predicts  that 
dynamic  temperature  variations  at  cold  and  hot  ends  under  zero 
cooling  load  agree  well  with  experimental  results. 

Three  types  of  operating  conditions  including  small,  medium, 
and  large  applied  currents  are  adopted  to  examine  the  effect  of 
material  properties  on  the  dynamic  characteristics  of  the  TEC. 
The  results  show  that  at  small  currents  the  constant  and  variable 
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property  models  predict  almost  the  same  dynamic  characteristics; 
at  medium  currents  the  predictions  of  two  models  differ,  and  the 
difference  is  more  apparent  at  high  cooling  load;  at  large  currents 
the  variable  property  model  predict  that  the  temperatures  at  cold 
and  hot  ends  increase  continuously  and  never  reach  the  steady- 
state,  oppositely,  the  constant  property  model  predict  that  the 
TEC  still  supply  refrigeration.  The  present  simulations  indicate  that 
the  variable  property  model  is  more  appropriate  to  accurate  pre¬ 
dict  TEC  transient  performance,  especially  when  the  TEC  operates 
at  medium  and  high  applied  currents. 

At  a  specific  applied  current,  the  temperatures  at  the  cold  and 
hot  ends  both  increase  with  the  elevated  cooling  load,  but  the  tem¬ 
perature  rising  rate  of  the  cold  end  is  much  higher  than  that  of  the 
hot  end,  resulting  in  a  significant  reduction  in  the  temperature  dif¬ 
ference  between  the  cold  and  hot  ends.  This  phenomenon  can  be 
attributed  to  that  small  temperature  difference  reduces  Fourier’s 
back  heat  transfer  and  hence  maintains  the  high  cooling  capacity 
of  TECs.  The  response  time  of  the  TEC  varies  from  several  hundreds 
of  seconds  to  several  thousands  of  seconds,  which  is  dependent  on 
its  operating  conditions. 
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